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Abstract 

A multifractal analysis is performed on the universality classes of random 
matrices and the transition ones. Our results indicate that the eigenvector 
probabi;ity distribution is a linear sum of two x^-distribution throughout the 
transition between the universality ensembles of random matrix theory and 
Poisson. 
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The Anderson localization is a wave phenomenon characterized by a destructive in- 
terference that gives rise to unaccessible regions in the configuration space of a physical 
system. It occurs in many situations [0 and, in particular, in condensed matter physics, 
it is the responsible for the metal-insulator phase transition (MIT) caused by an increas- 
ing of disorder in a quantum disordered system, in which, as a consequence, the material 
undergoes a transformation from a metallic to an insulator phase. The eigenstates of the 
system, in the metallic phase, extend over all the available space and, in the other hand, 
at the insulator side, it is localized around the impurities. This situation implies, as the 
transformation proceeds, in a modification of the fractal dimension of the wavef unctions. 
This aspect of the transition has been studied throw the use of the multifractal analysis |^,^ 
which was introduced some years ago p. 

Random matrix ensembles are another powerful theoretical tool to study transitions 
from extended to localized states 0. The phase of the extended states, i.e., the metallic 
phase in the MIT case, is approached by the universal ensembles of random matrix theory 
(RMT) 0], namely, the Gaussian Orthogonal Ensemble (GOE), if the system has time- 
reversal invariance, and the Gaussian Unitary Ensemble (GUE), if not. On the other hand, 
the insulator phase, where the states become localized, can be simulated by a Poissonian 
ensemble. Accordingly, the energy levels of the strong mixing metallic states are expected to 
follow Wigner-Dyson statistics of RMT, while the levels of the uncorrelated localized states 
of the insulator phase, have fiuctuations that follow the Poisson statistics. The Maximum 
Entropy Principle (MEP) has been used to generate matrix ensembles [|| that make the 
connection between two universal limiting situations, e.g., RMT and Poisson. The joint 
probability distribution of the matrix elements of the interpolating ensembles obtained are 
given by || 

P{H,f3,a) = Kr,exp[-aoTr{H')] exp{-(3Tr[{H - Ho)']} (1) 

where N is the dimension of the matrix, K]\f is the normalization constant, is fixed by a 
choice of units and /3 is the parameter that controls the transition. When (3 varies from zero 
to infinity, the ensemble undergoes a transition from the RMT ensemble to the ensemble of 
Hq. By choosing Hq to be a diagonal, we have the desired transition from RMT to Poisson. 
This connection between RMT and Poisson is not expected to be universal and, in fact. 



other possible interpolating ensembles have already been proposed, e.g., band matrices ||T0 
and U(N) invariant ensembles |]TT|. However, the above formalism has the advantage that 
the chaoticity parameter f3 is easily expressed in terms of the coupling constant In 
this letter, we extend the multifractal analysis to the abstract space of random matrices to 
investigate statistical properties of the eigenstates of this interpolating ensemble. 

Multifractal is a mathematical object that is not characterized by a unique dimension but 
by an infinite spectrum of dimensions [|13|. For example, in a classical dynamical system, 
the probability given by the frequency with which different "cells " of a partition of the 
strange atractor is visited, in the time evolution of a chaotic system, is a multifractal. In 
the quantum case, chaoticity may be defined as the situation in which the eigenfunctions 
spread uniformly over all components with respect to any basis. Of course, this refiects the 
fact that the system does not have any other good quantum numbers but the energy and it 
is the quantum mechanical equivalent of the absence of integrals of motion in the classical 
case. This property may be considered as an statement about the dimension of the wave 
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function. On the other hand, when transitions towards regularity occur, one should expect 
that the presence of conserved quantities implies in some localization that should be followed 
by modifications in the respective dimensions of the eigenstates. This can be seen in Fig. 
1, where it is plotted the logarithm of the components of a random eigenstate of the above 
ensemble as a function of its label, calculated at a critical value of the parameter /3 0. It 
is also shown in the figure, a Gaussian fit that makes clear the localization. It can be seen 
that only a small fraction of the components contribute to the normalization. In this sense 
the state does not occupy all the available space. Here the support is given by the states 
of the basis and the probability distribution by the way these basis states are populated by 
the eigenstates of the Hamiltonian. 

The critical behavior of the above transition ensemble has been defined in Ref. [^applying 
Shannon's concept of entropy to the eigenstates components treated as probabilities. As 
a start point to introduce the multifractal formalism, we observe that this entropy is a 
particular case of Tsallis generalized entropy WM of a probability distribution pi, with i = 



1, Np, which in terms of the partition function 

Np 

x{q)=T.Pl (2) 

i=l 

is defined by 

S, = i^. (3) 

for any real q. Associated with Sq, a spectrum of dimension functions, Dg, can then be 
introduced as 

0,^,^^ZS1=^, (4) 

^ «->0 (g - 1) In/ 

where / is a characteristic size associated with the partition. 

Some positive integer values of q have an immediate interpretation. Thus, Dq — 



Inl ' 

where is the number of occupied cells, i.e., those with probability different from zero, 
gives the fractal dimension of the support. For Di we obtain 

where 

Si = -'^Pi\npi (6) 

i 

is, by definition the Shannon's information entropy, and Di is the information dimension. 
D2 is the correlation dimension. 

We now assume the scaling pi ~ /" and, also, that the exponents a' vary continuously, 
as i is varied, with a certain distribution p{a') that scales with /■^^"^. Thus the partition 
function can be transformed into an integral as 
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X{q) = / da'p («') (7) 



Since / is a small quantity, this integral will be determined asymptotically by the value 
a' = a that minimizes the exponent of / in the integrand. This yields 

D, = [qa - f («)] (8) 

from which we can deduce 

« = ^ [(? - 1) D,] (9) 

These two equations give a and / (a) in terms of Dg or, alternatively, Dg if we know a 
and / (a) . Some universal properties follow from these definitions. Thus / (a) is a convex 
function whose maximum is located at a (0) with / [a (0)] = Dq and a varies in the interval 
[Doo, -D-oo] • / («) gives the dimension of the support of those pi that scale with a. 

As a first application of the above formalism, we consider the special case of the transition 
from GOE to the Poisson for matrices of dimension = 2. In this case, the probability 
distribution of a given component C can be worked out analytically. It has been shown in 



Ref. [15 1 that it is given by 




where < y = < 1. When j3 = 0, the GOE limit, the distribution is that of the 
component of a two dimensional unit vector that can point evenly in any direction. In 
the other limit, /3 — oo, the distribution goes to a sum of delta functions and the vector is 
completely localized, y can then have only the values or 1. This means that for small values 
of ratio (3/ao, the distribution is dominated by the two power law singularities located at the 



extremities of the segment. On the other hand, as /3 increases, the two poles M 1 =F . /l + ^ 
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of the denominator, approach the interval from the left(— ) and the right(+), respectively, and 
deform the power law behavior. In order to perform a more detailed analysis, we partition 
the interval [0,1] into Np equal subintervals of size I = 1/Np. We have then two kind of 
contributions to the partition function. Those that come from the power law singularities 
at ?/ = and y = 1 and those from the rest of the segment. To calculate the contributions 
of the first ones we integrate from to I and from 1 — / to 1 and for the others we just 
approximate them as p{y)l- The partition function is then given by 



xiq) 
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+ (11) 



For small values of the parameter /3, we can assume /3Z/ao -C 1 and the arctan can then 
be replaced by its argument and the first term becomes Z"^/^. We then deduce 

l2(g-l)' y ^ ^ 



4 



from which we derive a /(a)— spectrum with only two points, a = 1 with / (a) = 1 and 
a = I with / (a) = 0. These values mean that the extremity points have fractal dimension 
zero while the others have the dimension of the support. Incidentally, we remark that we 
have exactly here the same distribution of that provide by the iteration of the logistic map. 

In the second situation, we assume that although the / are very small (3 is so big that we 
can not linearize the arctan anymore. The first term scales then as P and we deduce 

^. = 1'' (13) 
\0, g>l ^ ^ 

and the two points a = 1 with / (a) = 1 and a = with / (a) = for the /(a)— spectrum. 
We remind that we are here getting close to the situation in which the distribution be- 
comes a sum of delta functions concentrated in the limits of the interval. The components 
approximate then the two values 1 and 0. That is why the exponent a vanishes. 

We see, in both cases discussed above, that the dimension function exhibits discontinuity, 
in its first derivative in one case and in the function itself in the other. In the thermodynamics 
picture of the multifractal analysis, these discontinuity are interpreted as phase transitions 
|16| and what we learn from the above results is that the transition chaos-order, here obtained 



by the variation of the parameter (3, is followed by a change in the qualitative behavior of 
the dimension function. This modification may also be considered as a phase transition with 
respect to the variation of the external parameter (3. This is exhibited in the figure where 
we can see that the information dimension Di plotted as a function of In ^ shows a typical 
first order phase transition pattern. One should expect that what we are observing in this 
case of ensemble of 2 x 2 matrices are universal features of the transition RMT- Poisson. 
We pass now to the discussion of ensembles of matrices of size arbitrarily large. 

We start discussing the GOE limit. It is known that, in this case, the probability 
distribution of the components is that of the components of a unit vector in the hypersphere 
in the space of N dimensions. It can be proved [IT^ that this is given by 



2r ( ^ ) 1 jv-3 
^fa)- „.,N„\i,N i'''(i-y)~ (14) 



The distribution is again characterized by the power singularities at the extremities y = 
and y = 1. The partition function associated to a division of the interval into Np cells of 
equal sizes / = 1/Np approaches, for small /, the behavior 

x(g) + /iv^. (15) 

and, in the limit when Z — > 0, we find the dimension function 

2 y ^ N-3 

D, = { 1, -^<g<2 (16) 

We have therefore a double phase transition separating three states, two of them are 
defined by the equation of state ^^Dq = const, which is generated by power law singularities 
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P and, in the middle of them, there is the state whose equation is given by Dq = const.. 
For the / (Q;)-spectrum we deduce, making use of Eqs. (|Dand @ , 

{0, a = I = -Doo 
1, g = l . (17) 

0, a = ^ = D^oo 

which means that the two singular contributions from the extremities of the interval have 
zero fractal dimension while the rest of the others have the dimension of the support. 

The more general case, with (3^0 and N > 2, has been investigated by numerical simu- 
lation of ensemble of matrices. The dimension function in the limit / was extrapolated 
from the dimensions obtained with two small I's by the relation 

_ D, {h)lnh-D, {h)\nh 
ln/i-ln/2 

which follows from the assumption that for sufficiently small /, the partition function behaves 

as 

The results obtained in this way for matrices of dimension = 100 are shown in Fig. 2. We 
see that the structure with two phase transitions and three states of the GOE case evolves 
to the picture given by Eq. (|TB|), typical of the limiting distribution with two 5-functions at 
the extrema of the interval. The correspondent / (a)-spectrum is shown in the next figure. 
Fig. 3, for three values of the chaoticity parameter and matrices of dimension N = 10. It 
is seen that the singular behavior of the GOE limit, Eq. ([T7|) has been smoothed out by 
the numerical simulation. Finally, in the subsequent figure. Fig. 4, the first-order phase 
transition exhibited by the information dimension Di, is seen. In the same figure, it is also 

2 

shown the Shannon's entropy of the eigenstates given by Eq. @ with pi = Cf averaged 
over the k states. In Ref. 0, the infiection point of this entropy has been taken as a definition 
of the critical value of the chaoticity parameter which separates the phase of localized and 
extend states. We can conclude from this figure that this definition is consistent with the 
results obtained for the information dimension. 

One important point to remark here is the consequence the present analysis has on the 
question of what is the probability distribution of wavefunction components and, also, of 
matrix elements of an operator- strength function- in the intermediate regime between RMT 
and Poisson or, in more general terms, between chaos and order. It has been proposed that 



a X distribution of v degrees freedom would fit this distribution |T8|. However, numerical 
simulations seem to suggest that a combination of two distributions is necessary in 



order to have a good description throw all the intermediate steps of the transition |]r2| . The 



above results of the multifractal analysis seem to give a theoretical support to this empirical 
observation. Indeed, what we have shown is that the chaoticity parameter acts like an 
external thermodynamic variable that induces a first-order phase transition.. Therefore, one 
should expect a modification of the nature of the probability distribution as the transition 
proceeds. In the N = 2 case, we have seen that this change in the structure of the function, 
Eq. (p!0D , is provided by the coming into action of two poles that lie outside of the physical 
domain. For large size of the matrices it is the appearance of an extra x^- distribution that 
takes care of the modification of the distribution in the passage from chaos to order. 
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Figure Captions: 

Fig. 1 Logarithm of the components of a random eigenstate for a case intermediate 
getween GOE and Poisson. The calculations were done with matrices of dimension N — 100. 
A Gaussian fit is also shown. 

Fig. 2 The dimension function Dg for matrices of size N — 100 for four values of the 

parameter P{ao = 1). 

Fig. 3 The / (a)-spectrum for matrices of size = 10 for four values of the parameter 
P{ao = 1). 

Fig. 4 Information dimension Di for matrices of size = 100, showing a first-order 
phase transition as a function of the logarithm of the chaoticity parameter. Also shown is 
Shannon's entropy of the eigenstates. 
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